Competing crystal structures in Lao.sCao.sMnOs: conventional charge order versus 

Zener polarons 
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Equilibrium crystal structures for Lao.sCao.sMnOa have been calculated using hybrid exact ex- 
change and density functional methods. Two distinct ground states with either conventional checker- 
board charge ordering or Zener polaron formation are found depending on the proportion of exact 



and d 3 y2 



exchange used. The checkerboard state has mixed d x 2_ y 2, d 3x 2 
CE-type magnetic order while the Zener polaron state has A Zx 2_ r 2 and d 3y 2_ 
magnetic order. 
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INTRODUCTION 

Low temperature charge ordering (CO) transitions in 
manganites such as Lai^^Ca^MnOs aregeneric for dop- 
ing in the range 0.4 < x < 0.9 0,11 EH 1 @. 
They are observed as commensurate or incommensurate 
changes in unit cell dimension parallel to the crystallo- 
graphic b axis |9j. Commensurate structures are found 
when x is a rational fraction such as 1/2, 2/3 or 3/4. Half- 
doped manganites, Ao.sBo.sMnOa with A = La, Nd, Pr 
and B = Ca, Sr, have been widely studied as the CO 
phase is believed to consist of a checkerboard (CB) pat- 
tern of Mn 3 + and Mn 4+ ions 0, E3, E3 in which Mn 3+ 
ions are Jahn- Teller (JT) distorted while Mn 4+ ions are 
not. However, there is evidence from Hartree-Fock calcu- 
lations [l2j, |l3( on Lao.5Cao.5Mn03 for a Zener polaron 
(ZP) electronic structure in which all Mn ions have a va- 
lence of 3.5 and a recent single-crystal neutron scattering 
study found a ZP crystal structure for Pro.6oCao.4oMn03 
[lij where all Mn ions have an intermediate JT distor- 
tion. There is conflicting experimental evidence from res- 
onant XRD experiments [15j for a CB CO pattern in 
P r o.6oCao.4oMn03 and a recent refinement of the struc- 
ture of Pro.sCao.sMnOs using high resolution x-ray and 
neutron powder data favored a CB structure ^3- Both 
CB and ZP phases have been proposed to exist in a phase 
diagram for half-doped manganites in which the tolerance 
ratio for the rare earth and alkaline earth ions is varied 
17|. Powder neutron and x-ray diffraction experiments 
and x-ray absorption near edge structure (XANES) 
flf| and resonant x-ray diffraction (RXD) [^3 sup- 
port a CB CO picture for Lao.sCao.sMnOs. 

In this Letter we use hybrid Hartree-Fock/density 
functional theory ab initio energy minimization calcu- 
lations and show that the ground state CO crystal struc- 
ture switches between CB and ZP as the percentage 
of Hartree-Fock (exact) exchange is varied. These two 
states are therefore similar in energy and the actual 
ground state may depend on specific A and B ions 0] , or 
even sample preparation, history and ambient conditions. 
The predominant difference in charge populations of the 



)",01a O 02 
_,},_Mn1a 



1*04 ; «03a #Qia ^ j> 



(a) 



it m 

|__lMn1_ 



^Mn2 



>04a ^03 



J 



I 



(b) 



01b Qo4 
A 03a 



(^Mn2 



(|)Mn2 I 

4 



4* 



T 



J 



FIG. 1: ab plane atomic positions for Lao.sCao.sMnOs shown 
to scale, (a) Pn2\m structure. Mn ion positions in the ZP 
state (Table are indicated by small circles, O 2- ions by 
shaded circles and 0~ ions in the ZP state by a filled circle. 
Mn and O ion positions in the checkerboard state (TableHJ are 
indicated by small and large dashed circles, (b) P2\nm ZP 
structure (Table ITU. Directions of transverse displacements 
of Mn ions from their ideal perovskite positions are indicated 
by arrows. Labels on ions refer to Tables U to Mil 



two states is in charge on O ions rather than Mn ions and 
so CO may be associated with O ions rather than Mn 
ions. Furthermore, we find that the lowest energy struc- 
ture for either CO pattern has a different space group 
from those proposed for the CB or ZP Q structures. 
The structure that we find has modulations of Mn and O 
ion positions about ideal perovskite positions which are 
parallel to ab planes and transverse to the b axis and 
are similar to those observed in layered [lllland simple 
doped manganites 0S0 00 A TA phonon soften- 
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ing mechanism for the CO transition would explain the 
observation of incommensurate order with wavevector q 
= (l-x)b* |(| and may also be relevant to colossal magne- 
toresistance (CMR) phases of manganites where diffuse 
satellites in x-ray scattering have been observed 0, 
and attributed to polarons with transverse displacements 
of ions (2^ similar to those in Fig. [Ja). 

Structure optimizations were performed using the 
CRYSTAL program ^ . A ^2x2^x2 unit cell (Fig. [TJ 
containing eight formula units with ferromagnetic (FM) 
order was used for optimizations while total energies of A 
and CE-type magnetically ordered structures were com- 
pared using 2v / 2x2v / 2x2 unit cells. A single force evalu- 
ation was carried out for CE-type magnetic order for a 
structure which had been optimized with FM magnetic 
order in order to estimate the importance of magnetoe- 
lastic effects on the crystal structure; these were found 
to be small. No simultaneous relaxation of unit cell di- 
mensions was performed as the calculations are expen- 
sive in computer time; unit cell dimensions were taken 
from experiment (a = 5.4763 A, b = 10.8932 A, c 
= 7.5247 A). Structure optimizations were performed 
for 60%, 80% and 100% exact exchange (the latter is 
simply an unrestricted Hartree-Fock (UHF) calculation); 
both 60% and 80% exact exchange calculations resulted 
in CB CO for all space groups investigated while 100% 
exact exchange resulted in ZP CO structures. Initial 
atomic configurations were taken from Table II in Ref. 
(P\Yl\l m symmetry) or were generated by hand. As 
noted by Daoud-Aladine and coworkers 01 > the isotropy 
subgroups of the parent high temperature Pbnm phase 
that have a doubled unit cell along the b axis are: P 11 m, 
P2 1 nm 1 P112 1 /m, P112i/6, P2 in b and P116. The rele- 
vant isotropy subgroups for ZP CO are P2\nm and PWm 
and for CB CO they are P112i/m and Pllra. 

Hessian matrices for energy minimized CB or ZP CO 
structures with these space groups had at least one neg- 
ative eigenvalue, which indicates that they are saddle 
points on the potential energy surface. A further energy 
minimization was performed using PY2\\ symmetry and 
100% exact exchange. The screw axis parallel to the b 
axis naturally incorporates transverse Mn displacements 
found in experiment. All Hessian matrix eigenvalues were 
positive for this ZP CO structure and it was lower in en- 
ergy than the other structures (Pllra +7 meV/Mn ion; 
P2 inm +14 meV/Mn ion). A structure optimization 
with 60% exact exchange and Pn2\m symmetry resulted 
in a stable CB CO structure which was 31 meV/Mn ion 
lower than the optimized structure with P112i/r7i sym- 
metry. Fractional coordinates for the lowest energy CB 
and ZP CO structures found, which have Pn2\m symme- 
try, are given in Table [I] Fractional coordinates for the 
Lao.sCao.sMnOa ZP CO structure are compared to those 
determined by neutron scattering for Pro.60Cao.40Mn.O3 
Q in Table El 

Fig. n^a) shows atomic positions for nearly coplanar 



TABLE I: Fractional coordinates for Lao.sCao.sMnOa with 
Pn2\m symmetry from ab initio total energy minimization. 
Wyckoff positions are given in the second column. Coordi- 
nates on the left were determined using 100% exact exchange 
and those on the right using 60% exact exchange. Both struc- 
tures have the (x,y) coordinates of the Mnl ion in common. 
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TABLE II: Fractional coordinates for Lao.sCao.sMnOa 
(LCMO) from ab initio energy minimization with 100% ex- 
act exchange and for Pro.6oCao.4oMn03 (PCMO) from single 
crystal neutron diffraction [Till , both with P2\nm symmetry. 
Wyckoff positions are given in the second column. 
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0.9021 


0.0000 
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0.0000 


La2 
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0.4783 


0.3723 


0.5000 
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0.3614 


0.5000 


La3 


2a 


0.9832 


0.1286 


0.0000 
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0.1426 


0.0000 


Ca4 


2a 


0.9741 
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0.5000 
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0.6088 


0.5000 


Mnl 


4b 


0.0000 


0.8743 


0.2498 


0.0000 


0.8756 


0.2489 


Mn2 


4b 


0.9611 


0.3775 


0.7518 


0.9795 


0.3746 


0.7492 


01 
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0.2717 
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0.2768 
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0.2861 
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0.2328 
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0.2276 
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0.4344 
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0.0000 


0.4353 


0.1125 


0.0000 
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0.5575 


0.1322 


0.5000 
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0.1321 


0.5000 
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0.0375 


0.3644 


0.0000 


0.0562 


0.3758 


0.0000 


04' 


2a 


0.9247 


0.3777 


0.5000 


0.9104 


0.3846 


0.5000 



Mn and O ions projected onto the ab plane in the Pn2\m 
energy minimized structure. Positions of ions in the ZP 
CO structure are indicated by shaded or filled circles and 
positions of ions in the CB CO structure are indicated 
by dashed lines (Table HJ). The structures are coincident 
at the Mnl ion positions and so this figure facilitates 
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0.0000 0.4861 0.8866 0.0000 

0.0000 0.4514 0.3673 0.0000 

0.5000 0.9743 0.1261 0.5000 

0.5000 0.9779 0.6057 0.5000 

0.2445 0.0000 0.8713 0.2478 

0.2534 0.4811 0.6231 0.2535 

0.2109 0.2352 0.9977 0.2109 

0.2112 0.2676 0.7600 0.2144 

0.2042 0.2118 0.2433 0.2146 

0.2753 0.7539 0.9970 0.2717 

0.5000 0.4320 0.1055 0.5000 

0.0000 0.5430 0.1337 0.0000 

0.5000 0.0316 0.3611 0.5000 

0.0000 0.9012 0.3625 0.0000 
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comparison of distortions in the ab plane. The main 
differences in ion position occur at the Mn2 and 03 posi- 
tions to accommodate the switch between ZP and CO JT 
distortion patterns. Transverse wavelike modulations of 
Mn and O ion positions about ideal perovskite positions 
can be identified in Fig. da). Modulations of either ion 
type have wavelength equivalent to the unit cell dimen- 
sion along the b axis and are out of phase by tt/2 and are 
similar to modulations of MnOfj proposed in manganites 
with x > 0.5 and x < 0.5 

Positions of ions in the energy minimized P2\nm struc- 
ture with ZP CO are shown in Fig. ^b). The main 
differences in atomic positions in the plane between this 
ZP CO structure and the Pn2\m ZP CO structure oc- 
cur at the 02 and 03 positions. Fractional coordinates 
for the P2\nm structures for both Lan.sCao.sMnOs (this 
work) and for Pr .6oCao.4oMn03 [lj| are compared in 
Table [H] and the agreement is remarkable. The starting 
guess for the Lao.sCao.sMnOs structure in ab initio cal- 
culations was not the Pro.6oCao.4oMn03 structure and 
the coordinates obtained via ab initio calculations are 
not simply a relaxation of the Pro.6oCao.4oMn03 struc- 
ture. There are bonds of intermediate length along the 
ZP axis ranging from 1.99 to 2.09 A in ab initio calcula- 
tions and from 1.98 to 2.05 A in experiment (Table 
llllll . Bond valence sums [24| (calculated using R Q and B 
values of 1.750 and 0.37 for all Mn-0 bonds) for the ZP 
structures are close to 3.5 in both computed and experi- 
mental ZP structures, which is indicative of an interme- 
diate valence. The Mn-0~-Mn bond angle is 160° or less 
in both Lao.sCao.sMnOs and Pro.6oCao.4oMn03. 




FIG. 2: Spin density and charge density difference for 
Lao.sCao.sMnOs in the ZP state with A-type magnetic order 
and Pn2\m symmetry obtained using 100% exact exchange, 
(left) Spin density. The unit cell shown in Fig. Q^a) is out- 
lined, (right) Charge density difference. 

The spin density and charge density difference of the 
ZP CO structure are shown in Fig. [2] Charge density 
difference plots are generated by subtracting densities of 
isolated 2 ~ and Mn 4+ ions from the total charge density 
of the crystal structures and therefore show deformations 
of charge density at O ion sites and e g orbital order at Mn 
ion sites. Each Mn ion in the ZP CO state has a d 4 con- 



TABLE III: Selected bond distances in A, bond angles for 
LCMO from ab initio energy minimization and for PCMO 
from single crystal neutron diffraction [Til ] and bond valence 
sums (BVS). Ion labels refer to atomic positions given in Fig. 
The bond angle given is the Mnla-03a-Mn2 angle at the 
center of the ZP or adjacent to a JT distorted Mn ion. The 
structure and space group are indicated at the top of each 
column. 



Bond 


CB a 


ZP 6 


Bond 


ZP C 


ZP d 


Mnla-Ola 


1.92 


1.99 


Mnla-Ola 


2.01 


2.05 


Mnla-04 


1.86 


1.87 


Mnla-Olb 


1.86 


1.88 


Mnla-02 


1.91 


1.89 


Mnla-02 


1.91 


1.91 


Mnla-03a 


1.84 


2.05 


Mnla-03a 


2.08 


1.98 


Mn2-02 


1.91 


1.89 


Mn2a-02 


1.90 


1.90 


Mn2-03a 


2.14 


2.09 


Mn2a-03a 


2.09 


2.01 


Mn2-01a 


1.92 


1.86 


Mn2a-04 


1.88 


1.90 


Mn2-04a 


2.15 


2.03 


Mn2a-04a 


2.02 


2.03 


Bond angle 


160 


155 




151 


159 


BVS Mnl 


3.28 


3.60 




3.60 


3.5 d 


BVS Mn2 
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a Lao.5Cao. 5 Mn03 Pn2 1 m symmetry, 60% exact exchange 
''Lao.sCao.sMnOs Pn2\m symmetry, 100% exact exchange 
c Lao.5Cao.5Mn03 P2\nm symmetry, 100% exact exchange 
d Pro.6oCa .4oMn03 P2 1 nm symmetry, expt. (Ref. [T^l 



figuration, O" ions order in Mn(d 4 )-0 _ -Mn(d 4 ) dimers 
(Zener polarons) in ab planes and the magnetic ground 
state is A-type 01 • UHF calculations 0] predict a mag- 
netic moment of 0.7 \xb on the O - ion in the center of 
each ZP, which is opposed to the moments of neighboring 
ZP Mn ions. On the other hand, cluster configuration in- 
teraction (CI) calculations [2j| show that the moment on 
0~ ions is much less than 0.7 fiB- However they show 
that the charge on these ions is approximately -l.Oe, in 
agreement with UHF calculations and that the ZP are 
strongly bound in a FM state. In both UHF and CI cal- 
culations the total magnetic moment on each ZP is 7/2 
Hb in agreement with experiment pj. The reason for the 
discrepancy between spin distributions in CI and UHF 
calculations is simply that the UHF spin function for the 
Mn-O - bond is a/3 whereas it should be (a(3 — (3a) /V2. 
d 3x 2_ r 2 and d3j,2_ r 2 orbital order in the ZP state can 
clearly be seen in Fig. [21 in both spin density and charge 
density difference plots. 

Bond lengths, the Mnla-03a-Mn2 bond angle and 
bond valence sums for the CB state with Pn2\m sym- 
metry are given on the left in Table II I II Spin densi- 
ties and charge density differences are given in Fig. [31 
The ground state magnetic order for the CB CO state 
is CE-type, the magnetic order found experimentally in 
Lao.sCao.sMnOs [2|, and it consists of zig-zag FM chains 
as shown in Fig. [21 Mnl ions in the chains have short 
(1.85A) Mn-0 bonds while JT distorted Mn2 ions have 
long Mn-0 bonds (2.15A). These are, respectively, the 
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FIG. 3: Spin density and charge density difference for 
Lao.sCao.sMnOs in the checkerboard state with CE-type mag- 
netic order and Pn2\m symmetry obtained using 60% exact 
exchange, (left) Spin density, (right) Charge density differ- 
ence. 

lengths of Mn-0 bonds in CaMn0 3 and LaMn0 3 USUI. 
The magnetic moments on these ions are 3.23/z# (Mnl) 
and 3.85/xb (Mn2) and bond valence sums are 3.28 and 
4.08. All of these are characteristic of conventional CO, 
although charges on Mn ions as measured by Mulliken 
populations are essentially identical, with values of +2.02 
(Mnl) and +2.02 (Mn2). A small difference in Mn ion 
charge has been noted in several experimental papers 
[lolllql in apparently CO states. 

In summary, we have shown that the equilibrium struc- 
ture for Lao.sCao.sMnOs which is predicted by hybrid 
Hartree-Fock/density functional theory depends on the 
percentage of exact exchange used in the calculation. The 
CB CO state consists of FM zig-zag chains in which cor- 
ner (Mnl) ions have short bonds (1.85A) while JT dis- 
torted (Mn2) ions have long bonds (2.15A) to neighbor- 
ing O ions. A preliminary calculation of the RXD spec- 
trum using this structure is in reasonably good agreement 
with experimental RXD spectra for Pro.6oCao.4oMn03 
\2l\ . The ZP CO state consists of FM planes tiled with 
polarons which have Mn-0 bonds of intermediate length 
along the polaron axis and this structure is in reasonable 
agreement with that found for Pro.6oCao.4oMn03 by neu- 
tron scattering [14j . Transformation from the CB to the 
ZP state requires rehybridization at Mnl and 03 sites to- 
gether with relatively minor displacements shown in Fig. 
^a). Effective magnetic moments from magnetic suscep- 
tibility data for at least two mixed valence manganites 
display marked increases in effective magnetic moment 
on cooling below the CO transition temperature to 6.1 
[lU or 7.9 hb These magnitudes are similar to the 
ZP magnetic moment of 7 /is . Dimerization of CE zig- 
zag chains to form a ZP state may produce this large 
effective moment. 
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